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USE  OF  THE  COMPUTER  PROGRAM  PARA  80  TO  STUDY 


RESULTS  FROM  FIRING  THE  RAILGUN  ERGS- 1 M 


1  •  INTRODUCTION 


This  paper  describes  simulation  calculations  of  railgun  performance 
using  the  computer  code  PARA  80  designed  by  Y.C.  Thio  [1].  The  calculations 
are  used  to  study  the  effectiveness  of  PARA  80  in  simulating  actual  railgun 
results  obtained  from  a  series  of  firings  of  the  ERGS- 1 M  railgun.  The 
program  was  originally  used  to  design  ERGS-1M,  so  the  performance  of  the 
latter  compared  to  the  predictions  is  of  much  interest. 


PARA  80  is  a  computer  program  written  in  FORTRAN  which  describes  the 
entire  railgun  circuit.  The  most  important  aspect  of  the  program  is  its 
description  of  the  plasma  armature.  This  is  done  with  a  one-dimensional 
model,  described  in  some  detail  in  ref.  1. 

The  hierarchy  of  subroutines  for  PARA  80  is  shown  in  Fig.  1 .  Table 
I  gives  a  brief  explanation  of  the  purpose  of  each  subroutine.  The 
calculation  may  be  considered  as  divided  into  two  parts.  In  the  first  part 
capacitor  discharge  dominates  and  in  the  second  part  the  capacitor  bank  has 
been  eliminated  from  the  circuit  and  inductance  dominates.  Consider  the 
railgun  circuit  diagram  shown  in  Fig.  2.  The  capacitor  bank  C  is  initially 
charged  to  a  nominated  voltage.  The  main  switch  SI  is  then  closed  and  the 
simulation  begins  at  this  point.  When  the  voltage  on  the  capacitor  bank 
drops  below  a  specified  value,  the  crow-bar  switch  S2  is  closed,  isolating  the 
capacitors  from  the  rest  of  the  circuit.  The  siaulation  then  switches  to  an 
inductively  driven  mode,  the  energy  of  the  circuit  being  taken  as  now  stored 
in  inductor  L.  The  program  ceases  calculation  when  the  circuit  current 
reduces  to  a  specified  value.  (Other  conditions  on  termination  and  crow-bar 
switching  are  allowed  by  the  program).  The  calculation  effectively 
progresses  down  the  barrel  from  the  plasma  Initiation  point  in  discrete 
(input)  timesteps.  The  distance  it  travels  down  the  barrel  is  determined  by 
the  projectile  motion  and  termination  condition,  and  may  vary  under  different 
simulation  conditions. 
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The  experimental  railgun  ERGS- 1M  ia  a  small  bore  (6x6  mm)  gun  with 
Cu-0.6%  Cd  rails  of  length  approximately  900  mm  [21.  The  insulation  between 
the  rails  is  a  vulcanised  cellulose  fibre  and  the  barrel  casing  is  made  of 
Micarta,  a  high-density  glass  reinforced  resin  material.  X  sketch  of  the  gun 
is  shown  in  Fig.  3,  full  specifications  are  given  in  Ref.  2. 


Instrumentation  of  firing  on  ERGS-1M  involved  recording  breech  and 
muzzle  voltages,  circuit  current  (with  a  Rogowski  belt)  and  projectile 
velocity,  in  the  barrel  using  inductance  coils  and  outside  the  barrel  using  a 
ballistic  pendulum. 

This  paper  reports  calculations  using  PXRX  80  to  compare  with  a 
series  of  firings  of  ERGS-1M  for  initial  capacitor  voltages  between  3.0  and 
7.0  kv.  The  simulation  program  was  run  on  a  VAX  11/780  computer;  Ref.  1 
described  the  program  structure  for  a  DEC  System  10  computer. 


2.  CHANGES  HEEDED  TO  RON  PARA  80  ON  VAX  1 1/780 


To  transfer  the  program  from  the  DEC  10  to  the  VAX  11/780  it  was 
necessary  to  make  several  changes.  In  order  to  make  the  data  file  containing 
the  run  specifications  more  legible,  it  was  reformatted.  This  was  done  for 
all  data  read  into  the  array  CON DAT  in  subroutine  SETUP;  the  nature  of  the 
input  data  is  shown  in  Table  II. 

In  addition,  the  subroutine  used  originally  to  solve  the  system  of 
equations,  NAG  library  routine  002ABF  has  been  removed  from  the  library.  The 
NAG  routines  on  the  VAX  11/780  are  in  double  precision,  and  subroutine  D02BAF 
was  found  to  be  a  suitable  replacement.  This  required  alteration  of  the 
arguments  of  the  subroutine,  particularly  accuracy  specification,  and  use  of 
double  precision  variables  where  relevant.  The  facility  for  reducing  the 
timestep  size  in  the  original  program  did  not  appear  to  be  very  effective,  and 
was  deleted.  A  check  on  the  timestep  size  is  now  done  through  changing  it  in 
the  input  data.  D02BAP  appears  in  subroutines  CAPSEG  and  INDSEG,  where  it 
uses  routines  CAP  AC  and  INDUCT,  respectively.  General  adjustments  to 
timestep  and  accuracy  parameters  were  also  required  to  operate  PARA  80  on  the 
VAX  computer. 


A  check  on  the  program  was  made  by  repeating  conditions  of  a 
calculation  previously  made  on  the  DEC  10  and  cosparing  results.  Although  it 
was  not  possible  to  reproduce  the  timeateps  exactly  as  before,  agreement  was 
sufficiently  close  to  accept  the  program  as  working  properly  on  the  VAX 
11/780.  The  disparity  arises  due  to  the  different  ways  in  which  the  original 
equation  solver  D02ABP  and  the  new  one,  D02BAF,  work.  The  former  would 
change  the  stepsize  if  it  needed  to,  while  the  latter  left  it  fixed.  This 
affected,  therefore,  the  positions  at  which  solutions  ware  found  along  the 
barrel.  As  a  result,  crowbarring  was  also  not  possible  at  exactly  the  same 
position  as  before. 
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3.  RESULTS 


The  experimental  firings  on  the  railgun  ERGS-1M  produced  records  of 
breech  and  muzzle  voltage  and  circuit  current  as  a  function  of  time.  Also 
recorded  were  projectile  positions  as  a  function  of  time,  and  from  these,  one 
can  estimate  velocity  as  a  function  of  time. 

PARA  80,  on  the  other  hand,  produced  curves  of  current,  position  and 
velocity  versus  time  directly.  Calculations  were  performed  at  the 
experimental  capacitor  voltages  of  3.0,  4.2,  5.0,  6.0,  7.0  kV.  For  each  of 
these  initial  conditions,  we  were  therefore  able  to  compare  current  profiles, 
but  more  interestingly  position  and  velocity  of  the  projectile  in  the  gun. 


The  calculations  were  repeated  twice  :  For  one  series  we  used  gun 
parameters  suggested  by  the  author  of  the  PARA  program.  These  are  shown  in 
Table  III.  For  the  second  series,  we  repeated  the  calculations  using 
parameters  thought  by  the  present  authors  to  be  more  realistic.  These  are 
also  shown  in  Table  III. 


To  effect  a  comparison  we  plot  projectile  displacement  versus  time 
for  each  of  the  firings.  Results  for  the  3.0  and  7  lev  cases  are  shown  in 
Figs.  4  and  5. 


The  positions  in  the  case  of  the  "experimental"  curve  are  determined 
from  magnetic  probe  measurements.  Ref.  3.  Six  probes  were  placed  on  top  of 
the  barrel  and  alternate  probe's  connected  together.  From  the  measured 
induced  voltages,  the  positions  versus  time  were  found.  The  accuracy  of 
these  measurements  is  subject  to  debate  since  it  was  found  that  the  voltage 
characteristic  of  adjacent  connected  probes  overlapped  to  some  extent.  This 
distorted  the  shape  of  the  characteristic,  and  added  to  the  difficulty  of 
interpretation  of  the  time  when  the  plasma  passed  the  probe.  Corrections  for 
the  difference  between  plasma  position  and  projectile  position  have  not  been 
made.  It  was  found,  however,  that  the  velocity  estimated  from  the  probe 
positions  agreed  reasonably  well  with  that  determined  from  the  ballistic 
pendulum  displacement  on  projectile  impact. 


The  results  of  Figs.  4  and  5  suggest  that  PARA  80  consistently 
under-predicts  the  projectile  velocity.  The  same  results  were  found  for  the 
4.2,  5.0  and  6.0  kV  simulations  as  well.  In  all  cases  projectile  velocity 
was  consistently  lower  than  experiment.  Changing  the  gun  and  plasma 
parameters  shown  in  Table  III  had  some  effect  -  generally  a  further  lowering 
of  projectile  velocity  -  but  did  not  dramatically  affect  this  overall 
conclusion. 
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4.  DISCUSSION  OF  THE  USE  AND  VALIDITY  OF  PARA  80 


PARA  80  has  been  used  by  Its  author  to  predict  the  performance  of 
various  firings  of  railguns;  the  predicted  data  on  armature  volt  drop, 
current,  and  projectile  position  as  functions  of  time  appeared  to  agree  well 
with  the  experimental  results  (1,4].  The  statement  is  made  in  the  "User 
Notes"  (1]  included  with  PARA  80  that  meaningful  interpretations  of  the 
results  obtained  by  running  the  code  requires  the  skills  of  a  plasma 
physicist.  It  is  well  known  that  computer  codes  tend  to  be  very  personal 
constructions  and  it  would  be  truer  to  say,  in  this  case  at  least,  that 
meaningful  interpretations  require  the  skills  of  the  code's  author.  Hie  User 
Notes  contain  such  statements  as  - 


PI  Has  the  value  of  1.0.  This  is  a  dummy  parameter  included  for 

compatibility  with  future  versions  of  PARA. 

A3  Dummy  parameter,  for  compatibility  with  future  versions  of 

PARA. 

A5,A6  These  should  have  the  value  0.45.  Included  for  compatibility 

with  future  versions  of  PARA. 

B6  Dummy  parameter 

B7  Should  have  the  value  of  1 .OE-9 

B8  Should  have  the  value  of  1 .OE+2 


Cl -Cl  1  [tammy  parameters  for  compatibility  with  future  versions  of 
PARA. 


AB  Should  have  the  value  of  0.0 

DU  Dummy  variable. 


Educated  guesses  can  be  made  about  the  purposes  of  some  of  these 
parameters,  but  the  author's  intention  about  others  is  not  clear.  In  the 
case  of  parameter  A6,  this  is  probably  the  same  as  y  (p.22  of  Ref.  1)  where 
it  is  stated  that  "1 -y_  is  the  fraction  of  current  lost  throuqh  arcing  ahead 
of  the  projectiles  normally  assume  y,  ■  1",  This  is  at  variance  with  the 
statement  in  the  User  Rotes  that  a  value  of  0.45  should  be  used.  In  any 
event,  there  is  no  physical  reason  to  assume  that  some  fixed  fraction  of  gun 
current  is  likely  to  pass  the  projectile. 

Re  believe  that  there  is  little  merit  in  combining  an  armature 
plasma  code  with  a  rail gun  performance  code  as  has  been  done  in  PARA  80, 
particularly  for  small  calibre  railguna.  The  performance  of  a  railgun  is 
dependent  on  one  plasma  parameter  only,  namely  the  resistive  voltaqe  drop 
across  it  (the  so-called  mutsle  voltage).  It  Bakes  more  sense  to  use  two 
codes,  one  for  railgun  performance  and  another  for  the  plasma.  It  is  known 
experimentally  (and  "confirmed"  by  computation)  that  mussle  voltage  is  largely 
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independent  of  current  for  a  particular  railgun  design  so  to  do  a  complete 
simulation  of  a  firing,  gun  system  performance  would  be  predicted  by  assuming 
a  constant  value  for  muzzle  voltage  to  give  current  as  a  function  of  time. 
Plasma  behavior  can  be  simulated  using  this  function,  to  give  in  turn  a 
computed  value  of  muzzle  voltage  as  a  function  of  time.  If  desired,  an 
iteration  could  then  be  made  for  gun  performance  using  this  computed  muzzle 
voltage  but  further  iterations  are  unlikely  to  be  required  because  of  the 
insensitivity  of  gun  performance  (including  gun  current)  to  muzzle  voltage. 

To  make  realistic  design  decisions  about  the  parameters  required  for 
future  railgun  systems,  simulation  codes  assuming  a  constant  value  for  muzzle 
voltage  are  adequate.  As  noted  above,  experimentally  it  is  found  that  muzzle 
voltage  is  essentially  independent  of  armature  current.  This  was  first 
observed  in  the  Canberra  railgun  (51  where  the  value  200  V  was  seen  with 
armature  current  ranging  from  300  kA  down  to  the  low  values  of  tens  of 
kiloamps.  The  same  constancy  was  observed  in  the  ERGS  experiments  [4]. 

Plasma  codes  may  then  be  used  to  enable  'fine  tuning”  to  be  done.  As  these 
simpler  codes  are  refined  by  comparison  with  experimental  results  it  should 
become  possible  to  construct  more  comprehensive  and  accurate  codes  which  could 
be  used  to  calculate  the  performance  of  large  calibre  railgun  systems.  Large 
systems  will  be  much  too  expensive  to  construct  without  having  very  advanced 
simulations  on  which  to  base  the  designs. 

5.  SIMULATED  PERFORMANCE  OF  AM  ERGS  FIRING  USING 
A  HP  65  PROGRAMMABLE  CALCULATOR 

To  demonstrate  the  adequacy  of  simple  simulations  for  small  calibre 
railguns  tie  include  the  following  model  which  can  be  programmed  onto  a  small 
calculator.  The  circuit  of  the  ERGS  system  is  as  follows  - 


L  R 


where* 


m 

v 

a 

MV 

x 

L' 


R' 


R 

L 

C 

V 

T 


I 
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is  the  railgun  projectile  mass, 
is  the  projectile  velocity 
is  the  projectile  acceleration, 

is  rail-to-rail  voltage  across  the  armature  (muzzle  volts), 
is  the  distance  the  projectile  armature  has  travelled, 
is  the  inductance  per  unit  length  of  the  railgun  rails, 
is  the  resistance  per  unit  length  of  the  railgun  rails, 
is  the  resistance  of  the  inductor  and  its  leads, 
is  the  inductance  of  the  inductor, 
is  the  capacitance  of  the  energy  storage  capacitor, 
is  the  voltage  produced  by  the  capacitor, 

is  the  voltage  drop  across  the  main  switch  when  it  is  conducting, 
similarly,  is  the  voltage  drop  across  the  crowbar  switch, 
is  the  current  flowing  in  the  railgun  circuit, 
is  the  time  rate  of  change  of  the  current. 


and  also 

t  is  the  elapsed  time, 

6t  is  the  time  increment  used  in  the  simulations. 


The  equations  (in  iterative  form)  which  describe  the  performance  of  the  system 
are 

-I  -  ( (R  ♦  R'x  ♦  L'v)I  ♦  (MV  +  T  -  V) )/(L  ♦  L’x)  (1) 

I  -  I  +  I$t  (2) 

a  -  L*I2/2m  (3) 


*  S.I.  units  are  used  throughout. 
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V 

*  V 

+ 

a6t 

(4) 

X 

*  X 

+ 

vit 

(5) 

t 

-  t 

+ 

6t 

(6) 

V 

■  V 

- 

I5t/C 

(7) 

After  the  capacitor  is  crowbarred,  (MV  +■  T  -  V)  in  equation  (1)  is  replaced  by 
(MV  +  T  +  S),  and  equation  (71  is  removed. 

Simulations  of  the  ERGS  shot  MIC9  (ref.  3)  have  been  made  with  the 
following  parameter  values 

C  »  0.00614  F 

L  *  6.8  E-6  H 

R  »  0.011  n 

MV  »  180  V 

T  «  S  »  2 OCV 

m  *  0.0003  kg 

R’  -  0.003  a 

V  initial  -  7000  V 

Crowbarring  was  assumed  to  occur  when  V  »  -500.  the  results  of  the  two 
simulations  with  I.*  _  0.25  E-6  and  L'  =  0.15  E-6  are  plotted  in  Fig.  5. 


6.  THE  CRITICAL  IMPORTANCE  OF  RAIL  INDUCTANCE 


It  can  be  seen  that  the  simulation  performed  with  the  lower  value  of 
L'  matches  the  experimental  data  well.  This  value  of  0.15  uH/m  is  very 
low.  The  value  for  L'  obtained  in  the  same  way  for  the  Canberra  railgun  (51 
is  0.42  uH/m. 

We  believe  that  the  difference  between  actual  performance  and 
measured  performance  has  much  to  do  with  frictional  drag  on  the  projectile  as 
it  is  being  accelerated  in  the  railgun.  Lower  values  of  L'  indicate  higher 
friction.  Higher  friction  effects  are  to  be  expected  in  smaller  guns  and 
this  is  what  the  results  show.  The  BIGS  railgun  has  a  bore  sise  of 
6  mm  x  8  am  compared  with  the  Canberra  railgun  bore  of  12.7  am  square. 
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It  can  be  argued  that  the  selection  of  some  particular  value  of  L' 
to  make  a  simulation  match  measured  performance  is  a  questionable  procedure. 
The  author  of  PARA  80  had  a  similar  factor  that  he  adjusted  to  get  agreement, 
namely  1  .  The  question  then  arises,  how  should  the  unknown  factors  be 
handled?2  The  answer  to  this  is,  "In  the  manner  which  gives  the  greatest 
insight".  Perhaps  then  we  should  be  simulating  by  using  a  computed  value  of 
L'  in  equation  (1)  and  modifying  the  acceleration  equation  (3)  to  read 


a  -  (0.5  L’  I2  -  F)/m 


where  F  is  the  friction  force.  If  friction  force  is  proportional  to  the 
driving  force,  then  the  effect  will  be  the  same  as  using  a  modified  value  for 
L'. 


It  should  be  noted  that,  in  any  case,  the  area  of  greatest 
uncertainty  in  the  simulation  of  railgun  systems  is  the  value  to  take  for  L' , 
the  rail  inductance  per  unit  length.  Depending  on  the  assumptions  made, 
values  can  vary  greatly.  The  three  main  possibilities  are: 


(1)  current  density  uniform  within  the  rails, 

(2)  current  confined  to  the  inner  surface  of  the  rails  to  the 
full  rail  height,  and 

(3)  current  confined  to  inner  rail  surface,  the  sheet  width 
being  equal  to  the  bore  height. 


For  the  Canberra  railgun  these  computed  values  [5,6,7]  are  0.503,  0.493  and 
0.628  uH/m.  For  ERGS,  the  respective  values  are  0.557,  0.552  and 
0.732  uH/m. 

Work  recently  reported  gives  hope  that  realistic  computed  values  for 
Ii'  will  be  available  soon  [8,9,10]. 


7.  CONCLUSIONS 


Fig.  5  contains  two  curves  estimated  by  using  a  simple  model  of  the 
railgun.  With  an  empirically  determined  value  of  the  rail  inductance  of 
L’  *0.15  v H/m,  good  agreement  is  found  with  experiment.  Dae  of  eqn  4.16  of 
ref.  1  indicates  a  value  of  L'  -0.8  pH/m,  which  is  considerably  larger.  In 
PARA  80,  it  is  thought  that  r.  of  eqn  4.16  [1]  is  input  rather  than 
calculated.  In  this  case,  a  value  of  t1  ■  0.45  as  recommended  in  the  "Oser 
Notes"  is  much  higher  than  the  value  of  0.24  obtained  from  Ref.  1  eqn  (4.17). 

Difficulties  in  interpreting  the  meaning  of  various  aspects  of 
PARA  80,  such  as  this,  considerably  limit  the  usefulness  of  the  program  for 


Jt 


the  present  work.  The  program  is  written  in  a  convoluted  manner,  and  uses 
unusual  scaling  of  variables,  both  of  which  make  it  very  difficult  to 
understand  from  the  source  listing.  The  possibilities  of  modifying  or 
improving  aspects  of  the  code  are  therefore  limited. 


We  believe  that  for  small  calibre  railguns  quite  simple  simulation 
codes  can  be  used  to  compute  performance.  Using  such  codes  and  experiments 
it  will  be  possible  to  determine  the  critical  physical  parameters  needed  to 
construct  more  comprehensive  codes.  Large  codes  will  be  needed  for  design 
and  analysis  of  large  calibre  railguns.  At  present  a  critical  area  in  need 
of  resolution  is  the  value  of  inductance  per  unit  length  of  rails  in  railguns 
and  this  is  being  addressed. 
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TABLE  I 


PAMPS  (MAIM) 
SETUP 

NORMAL 

EXTRAC 

SGUN 

INTERP 

CAPSEG 

DJDSEG 

RESULT 

J02BAF 

PHYCON 

MAGCON 

CAPAC 

INDUCT 

PLASMA 

VETROD 

VRIAL 


PURPOSES  OF  PARA  80  SUBROUTINES 


General  control,  files  opened  and  closed. 

Reads  input  data,  defines  constants  and  constructs  scaled 
quantities. 

Normalises  input  data  according  to  the  scales  from  SETUP. 

Structures  the  array  GUNDAT  for  the  actual  calculation. 

Controls  simulation  of  the  railgun  itself. 

Interprets  the  data  stored  in  GUNDAT. 

Controls  the  simulation  of  the  capactively  driven  gun. 

Controls  the  simulation  of  the  inductively  driven  gun. 

Prepares  data  for  output,  calculates  physical  properties 
and  outputs  results. 

Solves  (in  double  precision)  a  system  of  ordinary 
differential  equations. 

Calculates  pressure  behind  projectile  with  radiation  and 
magnetic  terms  only. 

As  for  PHYCON  but  includes  gas  pressure  term  as  well  (for 
closed  breech'  simulation). 

Calculates  the  integrand  for  D02BAF  for  the  capacitively 
driven  gun  -  i.e.  finds  the  array  F(10). 

As  CAPAC  but  for  the  inductively  driven  gun. 

Calculates  armature  plasma  properties. 

Gives  electrode  potential. 

Calculates  potential  drop  along  the  rails. 


-T* 
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TABLE  II 


-I 


I 

i 

THE  NATURE  OP  PARAMETERS  INPUT  TO  PARA  80  I 

| 


PARA  data  meanings 
Line  1 :  Type  of  foil 


Typical 

ALUMINIUM 


Line  2; 

Ionisation  potential  of  plasma  ions  (eV) 

Atomic  weight  of  plasma  ions 

Density  of  foil  Kg  m-3 

Height  of  foil  ) 

) 

Width  of  foil  )  metres 

) 

ThicKness  of  foil  ) 


5.98 
26.98 
2.7  E3 
6.0E-3 
44.0E-3 
1 5.0E-6 


Line 

_3s  CONDAT  (n,1) 

n 

i 

If  >  0.5  entirely  Inductive  simulation 

0.0 

2 

Reference  voltage  for  scaling  (Volts) 

1 0.0E+3 

3 

Reference  capacitance  (Farads) 

4.E-3 

4 

Another  dummy  parameter  (May  be  used  in 
routine? ) 

the 

NAG  file 

0.5E-6 

5 

Projectile  mass  (kg) 

.318- 

6 

Y,  -  0.45 

0.45 

7 

Y2  -  0.45  half  of  fraction  of  current 
through  arcing 

not 

lost 

0.45 

6 

Separation  between  the  rails  (a) 

8.8-3 

9 

Height  of  the  bore  (m) 

6.0E-3 

10 

Rail  height  (a) 

10.8-3 

i 

* 

i 

! 
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f 


Line  3  (contd ) 


Typical 


J 


11-20  Not  used  -  set  to  zero  0.0 

Line  4: 
n 

1  If  >  0.5  skips  norma lisation  (in  NORMAL)  (Unit  indicator)  0.0 

Indicates  whether  the  quantities  are  in  the  sane  line 

2  Electrode  potential  drop  per  electrode  (Volts)  10.0 

3  Magnetic  permeability  of  the  medium  (Hrn-1)  1.257E-6 

4  Resistivity  of  rail  material  (Qm)  2.E-8 

5  Height  of  the  bore  (m)  6.E-3 

6  Resistivity  of  armature  (Qm)  (Dummy  parameter)  4.E-8 

7  Width  of  the  bore  (m)  8.E-3 

_g 

8  Value  1.0  x  10  )  These  are  presumably  more  1.E-9 

) 

9  Value  1.0  x  10~*  )  dummy  parameters  0. 5E+2 

10  Resistivity  of  bus-bar  material  (Qm)  2.7E-8 

1 1  Length  to  width  ratio  of  bus-bar  1 60 

12  Resistance  of  the  circuit  other  than  bus -bar ,  rails  10.9E-3 

s  plasma  (Q) 

13-20  Not  used  -  set  to  zero  0.0 


Line  5:  CONDAT  (n,3) 

n 

1  ) 

2  ) 

3  ) 

4  )  Not  used  -  all  set  to  zero  0.0 

5  ) 

6  ) 

7  ) 
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Typical 


Line  6:  COHDAT  (n,4) 


n 

1  Number  of  segments  (should  be  1.0)  (NSEG)  1.0 

2  0.0  for  results  in  physical  (SI)  units,  1.0  for  0.0 

normalized  units 

3  Relative  error  of  computed  results  1.E-9 

4  No.  of  timesteps  skipped  before  computed  results  1.0 

output  to  files 


Line  7; 
n 

1  Unknown 

2  Unknown 

3  Inductance  of  storage  inductor  and  stray  inductance, 

excluding  rail  inductance  (H) 

4  Capacitance  of  capacitor  bank  (F) 

5  Initial  position  of  projectile  relative  to  closed 

breech  (si) 

6  Initial  current  (usually  zero)  (A) 

7  Initial  velocity  of  projectile  ( kins'"1 ) 

8  Initial  capacitor  voltage  (V) 

9  Initial  value  of  time  (usually  zero)  (sec) 

10  Tima  step  of  simulation  (sec)  (TINTLO,  TXNTVL) 

11  CONTYP  -  Controls  which  condition  is  used  for  crow-bar 

and  termination  (ICON) 

12  CONVAR  -  Criterion  to  control  crow-barring  and  final 

termination 


0.0 


o.o 


6.  6Z-6 


6.14E-3 


5.E-2 


0.0 


0.0 


3.E+3 


0.0 


1.3E-5 


8.0 


-500 
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Line  7:  (contd) 

13  CORSUP  -  termination  control  *  ^"^J^nt'hllff alien 

below  a  given  fraction  of 
the  current  at  the  contact 
of  crow-bar 


Typical 


14  Initial  Plasma  temperature  (K) 

15  ABSORB  (zero) ,  meaning  unclear 
16-20  Rot  used  -  set  to  zero 


2. 6E+3 


Line  8:  CORDAT  (n,6) 


Unknown 


Unknown 


3-15  Ro.  of  simulations  required  on  varying  the  relevant  (n,S) 
parameter 


16-20  Rot  used  -  set  to  zero 


Une  9:  CORDAT  (n.7) 


1  Unknown 

2  Unknown 

3.,5  of  <»;»  “ 

t(n,5)  gives  starting  value). 

16-20  Rot  used  -  set  to  zero 


not  used 


CHANGES  TO  INPUT  PARAMETERS  FOR  PARA  80  BETWEEN  THE 
ORIGINAL  AUTHOR  AND  THE  PRESENT  AUTHORS'  ESTIMATES 


CON DAT 

PARAMETER 

SIGNIFICANCE 

ORIGINAL 

NEW* 

(2,2) 

Electrode  potential  drop 

60  V 

10  V 

(4,2) 

Resistivity  of  rail  material 

4 

x  1 0“8  0m 

2  x  10“8  0m 

(10,2) 

Resistivity  of  busbar  material 

4 

x  10~8  flm 

2.7  x  10~8  Rid 

(12,2) 

Resistance  of  circuit  ex 
busbar,  rails,  plasma 

0.005  0 

0.0109  a 

*  Prepared  with  the  advice  of  Or  V.  Kovalenko 


- — Socfcat  Connection  _ . 

Low  Contact  Realatanca  AcWo*od 
By  Coopar  Wadgaa 


EROS-ttl  BARREL  AND  BREECH  ASSEMBLY 


FIG.  3 


Sketch  of  the  ERGS-1M  rallgun 


Plot*  of  position  va  tim  tor  initial  capacitor  voltage  of  3.02  kV. 
fa)  Rasults  of  the  author  of  PARA  80.  (b)  Results  for  the  "iaproved 
psraastsrs  shown  in  Table  III,  (c)  Results  as  for  (a)  except  for  a 
foil  asss  of  0.036  g*  (coapered  to  0,0104  ga  for  (a)  or  (b).  {d) 

Experimental  (MIC  11).  Ballistic  pendulua  MasursMnts  give  1500 
as"'  as  the  final  velocity.  Rough  estiMtes  fro*  the  curves  give 
the  following  velocities  near  the  end  of  the  barrels 


As  for  Fig.  4  but  at  an  initial  capacitor  voltage  of  7  kV 
Velocity  estimates  are: 

Ballistic  Pendulum  3300  ms  1  (MIC  2) 

(a)  187S  ms  ^ 

(b)  1450  ms" 

(c)  2050  ms  , 

(d)  3200  ms 

Curves  (e)  and  (f)  are  calculated  from  the  simple  model 
given  in  the  paper  for  (e)  L*  «  0.25 u  H/m, 

(f)  L'  -  0.15  y  H/m. 


-I 


(MRL-R-900) 


DISTRIBUTION  LIST 


MATERIALS  RESEARCH  LABORATORIES 


Director 

Superintendent,  Physical  Chemistry  Division 
Dr  A.  Bedford 
Dr  J.  Eadie 

Library  (2  copies) 

Dr  D.D.  Richardson 
Dr  R.A.  Marshall 


DEPARTMENT  OF  DEFENCE 

Chief  Defence  Scientist/Deputy  Chief  Defence  Scientist/  (1  copy) 
Controller,  Projects  and  Analytical  Studies/ 

Superintendent,  Science  and  Technology  Programme 
Army  Scientific  Adviser 
Air  Force  Scientific  Adviser 
Navy  Scientific  Adviser 

Of ficer-in-Charge,  Document  Exchange  Centre  (17  copies) 

Technical  Reports  Centre,  Defence  Central  library 
Central  Office,  Directorate  of  Quality  Assurance  -  Air  Force 
Deputy  Director,  Scientific  and  Technical  intelligence. 

Joint  Intelligence  Organisation 
Librarian,  Bridges  library 

librarian,  Bigineering  Development  Establishment 

Defence  Science  Representative,  (Summary  Sheets  Only) 

Australia  High  Commission,  London 

Counsellor,  Defence  Science,  Washington,  D.C.  (Summary  Sheets  Only) 
librarian,  (Through  Officer-in-Charge),  Materials  Testing 
Laboratories,  Alexandria,  NSW 
Senior  librarian.  Aeronautical  Research  laboratories 
Senior  Librarian,  Defence  Research  Oentre,  Salisbury,  SA 


DEPARTMENT  OF  DEFENCE  SUPPORT 
Deputy  Secretary,  DDS 

Head  of  Staff,  British  Defence  Research  c  Supply  Staff  (Aust. ) 


OTHER  FEDERAL  AMD  STATE  DEPARTMENTS  AND  INSTR (MENTALITIES 

NASA  Canberra  Office,  Woden,  ACT 

The  Chief  librarian.  Central  library,  CSXRO 

library,  Australian  Atonic  Energy  Commission  Research  Establishment 


MISCELLANEOUS  -  AUSTRALIA 


librarian,  State  library  of  NSW,  Sydney,  NSW 
University  of  Tasmania,  Morris  Miller  lib.,  Hobart,  Tas 


(MRL-R-900) 


DISTRIBUTION  LIST 
(Continued) 


MISCELLANEOUS 


Library  -  Exchange 
UK/CJSA/CAN/NZ 


Desk,  National  Bureau  of  Standards, 


(3  copies) 


Director,  Defence  Research  Centre,  Kuala  Lumpur,  Malaysia 
Exchange  Section,  British  Library,  UK 

Periodicals  Recording  Section,  Science  Reference  Library, 

British  Library,  UK 
Library,  Chemical  Abstracts  Service 

INSPEC:  Acquisition  Section,  Institute  of  Electrical  Qigineers,  UK 
Bigineering  Societies  Library,  USA 

Aero  medical  Library,  Brooks  Air  Force  Base,  Texas,  USA 
Ann  Germany  Documents  Librarian,  The  Centre  for  Research 
Libraries,  Chicago,  Ill. 

Defense  Attache,  Australian  Embassy,  Bangkok,  Thailand 
(Att.  D.  Pender) 


ADDITIONAL  DISTRIBUTION 


DARPA/TTO,  1400  Wilson  Boulevard,  Arlington,  VA  22209,  USA 
(Attn.  Dr  h.d.  Fair  Jr) 

ARRADOM,  Dover,  NJ  07801,  USA,  Building  382 
(Attn.  Dr  P.J.  Kammey,  O:  T.  Gora) 

ARRADCOM,  Ballistics  Research  laboratory  ERDAR-BLB,  Aberdeen 
Proving  Ground  MD  21005,  USA 
(Attn.  Dr  D.  Eecleshall,  Dr  J. D.  Powell) 

Naval  Surface  Weapons  Centre  (NSWC/DL) ,  Cahlgren,  VA  22448,  USA 
(Attn.  Dr  M.F.  Rose  Code  F-04,  Dr  H.D.  Odom  Code  F-12) 

Naval  Research  laboratory.  Code  4770,  Washington  DC  20375,  USA 
(Attn.  Dr  R.  Ford,  D:  I.  Vitkovitsky) 

US  Air  Force,  Armament  laboratory  (AFATL/DLCB) , 

Eglin  Air  Force  Base  FL  32542,  USA 

(Attn.  Dr  L.  Burge,  Dr  W.  Lucas,  Lt  R.  Walley) 

US  Air  Force,  Aero  Propulsion  laboratory.  Power  Division  (AFWAL/POOS-2) 
Wright-Patterson  Air  Force  Base  OR  45433,  USA 
(Attn.  Cr  C.  Oberly) 

US  Air  Force,  Office  of  Scientific  Research,  Bolling  APB, 

Washington  DC  20332,  USA  (Attn.  Cr  L.  Caveny) 

US  Air  Force,  Rocket  Propulsion  laboratory,  Edwards  AFB, 

CA  93523  (Attn.  LKDH/Lt  Phil  Roberts) 

NASA  Ames  Research  Centre,  Moffet  Field  CA  94035,  USA 
(Attn.  Technical  Library) 

NASA  langley  Research  Center,  Langley  Field  Station,  Hampton  VA  23365, 
USA  (Att.  Technical  library) 

NASA  lewis  Rssserch  Center,  Cleveland  OH  44135,  USA 
(Attn.  William  R.  Kerslake) 


(MRL-R-900) 


DISTRIBUTION  LIST 
(Continued) 


ADDITIONAL  DISTRIBUTION  (Contd) 

NASA  Goddard  Space  Flight  Centre,  Greenbelt  MD  20771,  USA 
(Attn.  Technical  Library) 

NASA  Scientific  and  Technical  information  Facility, 

Baltimore  MD  21240,  USA 
(Attn.  Accessioning  Department) 

Jet  Propulsion  laboratory,  4800  Oak  Grove  Drive,  Pasadena  CA  91102 
USA  (Attn.  Technical  library) 

Batelle  Columbus  laboratories,  505  King  Avenue,  Columbus  OHIO  43201 
USA  (Attn.  Dr  Eric  Rice) 

General  Dynamics,  Kearney  Mesa  Plant,  P.O.  Box  1128,  San  Diego 
CA  92112,  USA  (Attn.  Cr  J.H.  Cuadros  (Pomona  Division)) 
International  Applied  Physics,  Research,  7546  McEwen  Road, 

Dayton  OHIO  45459,  USA  (Attn.  Er  J.P.  Barber) 

Lawrence  Livermore  National  laboratories,  PO  Box  808, 

Livermore,  CA  94550,  USA  (Attn.  Cr  Ron  Hawke) 

Litton  industries,  360  North  Crescent  Drive,  Beverly  Hills  CA  90210 
USA  (Attn.  Dr  R.  Salter,  Dr  J.  Scudder) 

Los  Alamos  National  laboratory.  Mail  Stop  525,  Los  Alamos  NM  87545, 
USA  (Attn.  Dr.  Jerry  Parker) 

Los  Alamos  National  Laboratory,  Mail  Stop  985,  Los  Alamos  NM  87545, 
USA  (Attn.  Dennis  Peterson) 

Massachusetts  institute  of  Technology,  Francis  Bitter  National 
Magnet  Laboratory,  Cambrdige  MA  02139,  USA 
(Attn.  Cr  H.  Kolm) 

Physics  international  Company,  2700  Merced  Street,  San  Leandro 
CA  94577,  USA  (Attn.  CT  Ed  Goldman) 

Science  Applications  be.  6600  Powers  Ferry  Road  (Suite  220), 
Atlanta  GA  30339,  USA  (Attn.  Dr  J.H.  Batteh) 

University  of  Texas  at  Austin,  Center  for  Electromechanics  (CEM), 
Austin  TX  78712,  USA  (Attn.  Or  William  Weldon) 

Vought  Corporation,  Advanced  Technology  Canter,  PO  Box  226144, 
Dallas,  TX  75266,  USA  (Attn.  Or  C.H.  Haight,  Cr  Mike  Tower) 
Westinghouse  Electric  Corporation,  1310  Beulah  Road, 

Pittsburgh  PA  15235,  USA 

(Attn.  Mr  C.J.  Mole,  Cr  I.  Me  Nab,  cr  D.W.  Deis,  Or  Y.C.  Thiol 
Westinghouse  Electric  Corporation,  Sunnyvale  CA  94088,  USA 
(Attn.  Mr  William  A.  Vols) 

Sand  la  Laboratories,  Albuquerque,  NM  87115,  USA 
(Attn.  Cr  M.  Cowan) 

Rand  Corporation,  Santa  Monica,  CA,  USA 

(Attn.  Dr.  Gerald  A.  Sears,  Cr  William  J.  Whelan) 

General  Atomics,  San  Diego  CA,  USA  (Attn.  Cr  Sibley  Burnett) 

JAYCOR,  205  South  Whiting  Street,  Alexandria  VX  22304,  USA 
(Attn,  cr  Derek  Tidman) 

General  Electric  Company,  Schenectady  MY,  USA 
(Attn.  Mr  Robert  A.  Marshall) 

Research  School  of  Physical  Sciences,  Australian  National 
University,  Box  4  PO,  Canberra,  ACT  2601 
(Attn.  Prof.  S.  Kaneff) 


* 


END 


DATE 

FILMED 


